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1 Introduction. 

The problem of averaging for systems with one fast phase was considered from various points of 
view in many papers. The averaging method of Krylov and Bogolyubov |jl| and methods of KAM 
theory originated this line of research, the most complete results were obtained by Neishtadt [Q], 
where the coefficients are assumed real analytic. However, in many problems which are interesting 
from the point of view of applications, analytic dependence ceases to hold in the neighborhoods 
of some points. 

For instance, consider the motion of a particle under the influence of a small periodic electric 
and constant magnetic field [0]. The Hamiltonian of such a system has the form : 



The canonical change of variables xi = Q + yi, pi = — 1/2, X2 = P + 1/2, P2 = —Q, Q = v 2/ cos 99, 
P = a/27 sin transforms it into: 



H = I + eV ^-\/27cos(y9 + v^27siny9 + 7/2 j 



This system depends on one fast phase ^ and does not depend analytically on / near / = 0. This 
is not a minor problem since for many interesting problems coming from physics, a neighborhood 
of / = may play the most important role. In the problem of semiclassical quantization, for 
examples it corresponds to the so-called low levels of Landau which are connected with Hall's 
conductivity. This problem is actually the main motivation of present paper. 

The procedure in p is based upon subsequent change of variables, which corresponds to 
methods of KAM-theory. The aim of this paper is to show that one can choose a transformation 
such that the averaging procedure of [0] is applicable in a neighborhood of / = 0, and such that 
the passage from / >x>0to/ = 0is uniform. 



1 



2 Formulation of the problem and the main result. 

Consider the Hamiltonian 

H = no{I) +6go{q,p,yi,y2), 1= — - — , (1) 



where < e < Eq is a small parameter, and T-Cq and Qq are real analytic functions in a complex S- 
neighborhood of the domain D := D2n{yi,y2} x C M^", = {{q,p) G | / < x}. 

In D we assume the following conditions 



l^nl < C, 



\go\ < c, 



dl 



Following |1|, 0], we show that for every integer m > there exists a close to identity real 
analytic canonical transformation {q,p,yi,y2) ^ {Q, P, Zi, Z2) defined by 



where + |P^| 



q = Q + sQ^iQ, P, zi, Z2, e), y^ = zi + eZ^Q, P, zi, Z2, e), 
p = P + eP^{Q, P, zi, Z2, e), y2 = Z2 + eZ^iQ, P, zi, Z2, e), 

\Zl \ + 1^2! < C, and such that the Hamiltonian ( |l| ) transforms into 
'Q' + P' 



(2) 



H = nr 



, Zi, Z2,s + egm{Q, P, Zi, Z2, e) 



(3) 



Precisely, we have the following theorem. 

Theorem 1. Assume the conditions above in {Q, P, zi, Z2) G -D + ^S. Then there exists some 
interval (0, ei], integer number r and real analytic canonical change of variables of the form( |^ ), 
which transforms Hamiltonian ( |^ ) into ( ^ ) with exponentially small g: 



\gr\ + \^gr\ < C2 exp 



CiE 



\Q^\ + \Zl\ + \P^\ + \Zl\ < C3, \Hr - Hoi < C4£. (4) 



Here e G (0, ei] and 5i, r, q, z = 1, 2, 3, 4 depend on Eq, 6, C and x. 



3 Auxiliary lemmas. 



Let w (/, /i) and g{q,p, jj) be an analytic function of / and respectively and vector-parameter 

dq ~ ai' dp ~ t' dr 

dw dW dw dW 



^ ^ J = (g2 +p2)/2. Denote ^ = q^, ^ = p^. Consider the equation : 



dp dq dq dp 

_q^+p^ 
9{ — — ,/^) = 



g2 _^ p2 



2lT 



'0 

-,2 I ^2 



{l=I(q,p)) 



( q \ p 

9{<l,P,f^) = 9 — 7^ — ] -9{(l,P,f^), 



(5) 
(6) 

(7) 
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where g{{q'^ +p^)/2,/i) is a mean value of g{q,p,ij) with respect to ip, and the rest part of g is 
noted as g{q,p,fi) and it is convenient to take sign minus here. 

Lemma 1. Equation (H) is solvable and has analytic solution on variables (g, p) and parameter 
/i. Function W{q,p,fi), defined by formula : 



^ j 9{Q{i^J),p{'4'J),fJ')d^+^ j g{q{i),I),p{^,I),lJ,)dip 



is an analytic solution of equation (|^). 

Remark. General solution of equation IM) has a form : 



(8) 



W{q,p,fi) 



dw 
dl 



giq{ip, I),piip, I),ix)dip + Woil, /i) 



v>='Pi<i,p)J=Hg,p) 



where Wo{I,iJ,) is a constant of integration. It follows for {q,p) variables from method p|. We 
show below, that in the neighborhood of point q = 0,p = definition of constant of integration, 
pointed in ( H ), allow us to integrate equation ( ^ ), preserving analyticity of solutions. 



Proof. 



First make a canonical change of variables: 

u + iv 



q 



Then equation (13) is : 



V2 



dw dW dw dW 



P 



V + iu 
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u + iv V + iu 



dv du du dv ) ^ \ ' ^/2 
Now we should expand g into Taylor series. The expansion has form: 



(9) 



u + iv v + iu \ ^ , , u Ak + l)\ 

9\ ^,/^ ) = 2^ ^H(/i)MV' 



k\l\ 



Using coordinates {(p,p): u = pe*'^, v = pe one can easily integrate equation 



W 



dw 

dp 



u + iv V + iu , , , , , / N 
9 I — — p]dil) \ + PVolP, N 



^ 9 klip) 

dp \k,l€N 



V2 ' V2 
u^{^,p)v\^,p) (k + iy. 



i{k-i) km 

where Wo{p, p) is a constant of integration. 



9m{p)-7t- 



(k+l) 



{k + l)\ 



k,ien 



i{k-l) kill 



+ Wo{p,p), (10) 
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We can see that at bottom limit of integration we can have nonanalytic ( uv = v27 ) 
dependence on / at the point {u = 0,v = 0). But we can choose Wo{p, p) in a such way that 
nonanalytic term disappears. Introduce Wq as : 



dp \k,im 







l{k-l) kill 



Hence ([TT]) 'kills' nonanalytic term in (PUj), so W is analytic function of I in the point / = 0. It is 
easy to see that ^ = |j. Function W{q^p, p) is analytic ones of parameter p due to procedure of 
construction: we had no nonanalytic dependence on p on all steps of obtaining of solution. Hence, 
W is an analytic function of and parameter p. Lemma is proved. 



We use generating function to construct canonical transformation in the proof of main theorem. 

Let S{q, P,yi, Z2,e) be an analytic function of all variables {q, P,yi, Z2) and small parameter 
e e [0, £0], {q,P, yi, Z/2) e U, {Q, P, z^, Z2) e U - S, S > 0, U C M^n+z, Domain U - 5 is the set of 
points from U, which enter together with their 5- neighborhood |^ . Consider a system of equations 

dS{q,P,yi,Z2,e) dS{q, P,yi, Z2,e) 

Q = q + e — , zi = yi+e 



dP dz2 .^2) 

dS{q,P,yi,Z2,e) dS{q, P,yi, Z2,e) 

p = P + e , y2 = z2 + e , 

dq dyi 



f \dS\ 


1 1 


1 as 1 


1 1 


\ 1 1 ' 




1 dyi 1 





Following lemma establishes that introduced transformation is a change of variables and allows 
us to estimate corrections. 



Lemma 2. If e < 5/i2Cu{n + 1)), where n is a dimension of j/i, then the system ([12D is 
solvable and solution has the following form: 

q = Q + eQ^{Q,P,zi,Z2,e), yi = zi + eZKQ, P, zi, Z2,e), ^^^^ 
p = P + eP^{Q,P,zi,Z2,e), y2 = Z2 + eZ^iQ, P, zi, Z2,e), 

where , P^ , Z\, Z]^ are analytic functions in f/ — 5 x [0,60] and max{|(5^|, \P^\, \Z\\, 1^2 1} < 
in domain U — 6. 

Proof. 

Let us present system ( [T^ ) as : 

fiiq,P,yi,y2;Q,P,zuZ2,e) = 0, f3{q,p,yi,y2;Q,P,Zi,Z2,e) = 0, ^^^^ 
f2{q,P,yuy2;Q,P,zi,z2,e) = 0, f4{q,p,yi,y2;Q,P,zi,z2,e) = 0, 
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and F = (/i, /2, /s, /4). Due to the theorem of imphcit function, the system ( p^ is solvable with 
respect to {Q, P, Zi, Z2), if det I D{qpyi ^2) I ^ ^^"^ solution is analytic in f/ — 5 x [0, ^o]- In our 



case 



det 



DF 



D{q,p,yi,y2) 



det 



1 + e 



dPdq dPdyi 



8x28 q 



l + e 



dz2dyi 



S{q, P,yi, Z2,e) is an analytic function of e. Using Cauchy estimations for analytic functions [^], 
we obtain 



det 



DF 



D{q,p,yi,y2) 



71+1 



i=l 



s 



{n + iy. 

(n + 1 



< 



<1. 



(15) 



Hence, ii e < 6/ {2Cu{n + 1)) then ( |I5| ) is valid on U always and our system is solvable. 
Uniqueness of the transformation is established similarly to The solutions should satisfy 
the equations (|12D with e = also, so they can be represented as (|1^). Let us substitute the 
solutions in equations (|I2p. Then from equations for q and yi we have the following form for 
and Zl: 



dS{Q + eQ\P,zi+eZlz2,e) 



dS{Q + 6Q\P,zi + 6Zlz2,e) 



dP ' ' dz2 

Hence \Q^\ < Cu and | < Cu- For p and y2 variables we have the following formulae: 



p = P + e 



dS{Q + 6Q\P,zi+6Zl,Z2,e) 



y2 = Z2 + €- 



dq 

dSiQ + eQ\P,z, + eZlz2,e) 
dyi 



P + e 



Z2 + e 



dS{Q + eQ\P,zi+eZlz2,e) 
dQ 

dS{Q + eQ\P,zi + eZlz2,e) 
dzi 



Therefore: 



,1 ^ dS{Q + 6Q\P,zi+6Zlz2,e) 
dQ 



,1 _ dS{Q + 6Q\P,zi + 6Zl,Z2,e) 
dz. 



and again |-P^| < Cjj, {Z^l < Cu. Lemma is proved. 



4 Proof of the theorem. 

We construct our change of variables as a composition of large number of consequently defined 
canonical transformations, giving dependence on {q'^+p'^)/2 of Hamiltonian more and more higher 
degree of e. 
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4.1 Procedure of consequently defined changes of variables. 

Assume that Hamiltonian is obtained after i changes of variables has the following form 



H = Hi 



q2 _j_ p2 



yi,y2,£] + egi{q,p,yi,y2,e), 



(16) 



{(i,y2,p,yi) e A, 

D + -cD,cD + 5, 



Di = Di~ 2(i - l)Ke, 
Di = D + ^6. 



At i + 1 step one has to find canonical infinitesimal ( "almost identical" ) change of variable (see 
II, i) iq,p,y,,y2)^iQ,P,z,,Z2) 



q = Q + eQ^{Q, P, zi, Z2, s), y^ = zi + eZ\{Q, P, zi, Z2, e) 
p = P + eP^{Q, P, zi, Z2, e), y2 = Z2 + eZl{Q, P, zi, Z2, e) 



(17) 



such that Hamiltonian takes a form : 

'g2 + p2 



H = 'H 



zi,Z2,e) + egi+i{Q, P, zi, Z2, e) 



where 7ii+i contains terms of order i + 1 of e, and the operation "bar" over Qi is defined in sec. 3. 

It is possible and convenient to do by means of generating function (then one has the canonical 
property automatically) : 

S = S{q,P,yi,Z2,e) = qP + yiZ2 + eS^{q, P,yi, Z2, e). 

All other variables are defined in terms of (g, P, yi,Z2) and S{q, P,yi,Z2): 



Q = q + e 



dS^{q,P,yi,Z2,e) 
dP 



zi = yi+e 



dS\q,P,yi,Z2,e) 

p = P + e , y2 = Z2 + e 

oq 



dS^{q,P,yi,Z2,e) 
dz2 

dS^{q,P,yi,Z2,e) 
dyi 



(19) 



Hi 



Hi 



Substitute (|T^) to (|T^) and (|T8D, and equate Hamiltonians in "mixed" new-oldcoordinates : 



yi,Z2 + e——,e 
oyi 



\ 



\ 



,yi+ ^^^,Z2,€ 
OZ2 



+ egi{q,P + e—-, yi, Z2 + e——, e 
oq dyi 



q ^ t. gp 



J 



\ 



,yi + e——,Z2,e 

OZ2 



+ 



+ EQi+i \ q + e-^, P, 2/1 + e-^, Z2, e 
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Expand them in Taylor series, write out the terms of the same degree of e and take into 
account, that dHi/dyi = dT-Ci/dz2 = 0{e) : 



T-^i I — 2 — ' ^ ) =^i[ — 2 — ' ^ ) ' (20) 

+ 9i{Q, P, yi, s) = gi[ , yi, Z2,e ] . (21) 



dP dq dq dP 

Equation (^) is integrated with Lemma 1: S{q, P,yi, Z2,e) = W{q, P, fi), where n = {yi,Z2,e). 
Function S is analytic function of all its variables and small parameter e, so we can solve system 
([19|), using Lemma 2. The solution is continuous and has continuous derivative, so it defines 
change of variables. We find necessary canonical transformation (p!7D and after substitution into 
([T6|) we can obtain explicit expression for 7ij+i and egi^i. 



4.2 Estimations. 



Let r steps have been done. Domain Dj, where our Hamiltonian is considered after i steps, defined 
as A+i = Di — 2{i — V)K£, where = Z) + |5, K -positive constant, defined below. 

At first step we can easily show, that if {q,p,yi,y2) ^ Di and e is quite small, then ([17D is 
defined and the following conditions are realized: 



1^1 1 



\\/gi\ < he, 



+ P^ + Zl + Zl\ < k 



2, 



|7Yi — TYol < he. 



{ki] are positive constants. Really, from ( [T^ ) we immediately obtain that gi = (^(e^). Using 
estimation procedure defined below all other inequalities are easily derived. 

Let us accept inductive hypothesis that with i : 1 < i < r following estimations hold: 



m < 2c, 



c < 



dHi 



dl 



< 2C, 



\Vgi\ + \g^\< M,, 



Mi 



kie 



-1 ' 



(22) 



where / = (g^ +p'^)/2. Now one should find ei and K such that if < e < ei, {Q, P, zi, Z2) G 
Dr+i = Dr — 2Ke, the system ([T7|) is defined with i = r and ( p2l) is satisfied with i = r + 1. 

For we have the following estimations from definition of ( ) j equation ( ^ ) and form 
of solution ( Bl ) : 



OS' 



dip 



< \g,\ < \gi\ < Mi. 



The estimations for J = [q^ + P'^) /2 is obtained by the following way : 




QJ2 ' an^dj 

aJ aJ 



giip, J)di) + - / g{i), J)d'4) 



- / g{i), J)di) + - / g{i), J)di) 



(23) 
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We can estimate ^^jr- by : 



< 



+ e 



Using initial data we can state that l^jl^l < ""^i and \ ^jt\ < k^, where {rrii} are positive constant 
and doesn't depend on step i. For 



I aj2 1 



Summing up, we obtain 



< 



we use Cauchy estimations for analytic functions 



< 



EkiE 2ki 



9J2 



2A;i 

< mi + fc4£: + — £ < fcs. 



(24) 



For ^ we have the following estimations from (E2| 



dJ 



dp 



dJ 



dgi 



dp 



1 

< - 
P 



dgi 



dP 



+ 



dq 



< -2Mi. 
P 



where p^ = J. 

Then integral in ( ^3]) can be estimated as: 




dJ ^ 2 



d~g{iJ,J) 
dJ 



dip 



An Mi. msMi 
< ) < 



P 



P 



Then expression ( |23| ) is estimated: 



dJ 



< {k, + ^)M,. 
P 



Now we have the estimations for ^ and ^ and can obtain estimations for ^ and 







dP 






< 




< 



dJ dS^ dip dS^ 
+ 



dP dJ dP dip 



< p{k, + ^)Mi + 
P 



dgi ^ dgi p 
dq dP 



< 



m < 



(25) 



For estimations obtain the same way. 



dS^ 
dq 



< m^Mi. 



(26) 
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C-L 

Estimations for the ^f- and jf- are obtained easily by differentiation of Oi in the definition of 

ayi OZ2 J J 

solution for 5* : 



dyi 



< rriQMi, 



OS' 



< mjMi. 



(27) 



Now we can define fcg = niax(m4, ms, mg, rrij). Function gfj+i can be defined by : 



\gi+i\ 



2 \ dPy g \ dq J 2 \ dq^ J g\dP 



+ 



+e 



^g^\ dS^ fdg,\ dS^ f dgA dS^ 



dP J n dq 



— e 



dz2 J g dyi V 9yi J g dZ2 
^g^\ dS' 



— e 



dq J dP \dz2Jndy1 ydyij^dz. 



where {)g means the derivative in the middle point. Thus, gi+i estimates using (0), (p^)-(|27|) : 



\9i+i\ < 



+ 2e\ki + 



2ki 



Mi + AekaMt 



< kjM^.e 



Using Cauchy estimations we obtain estimations for Vgi+i. 



\^9i+i\ < 



4- 



ek^keM? 



2Ke 



+ 2e\kA + 



2ki\ Mi 



K Ke 



Ml 
'Ke 



< ksM, e + 



K 



Choosing K quite large and e quite small we obtain kje < \^ k^ [e + < \. Then 



, ,^ , Mi _ 

\9i+i\ + \ygi+i\ < ^ = Mi+i, 



and other inductive inequalities (p2|) are valid for i = r + 1. Therefore, we can do necessary 



changes of variables with chosen K,e till Dr is not empty. After r = (jl^) > kg/e changes we 



have obtained: 



II 1,-^ I / 1 N 

\gr\ + V^r < — - < C2exp{ ), 

2'^~^ CiS 



and other inequalities (^) will be satisfied. 
Theorem is proved. 
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